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Abstract 



■ In Ref. we presented an asymptotic formula for the fermion-induced effective 

energy in 3+1 dimensions in the presence of a cylindrically symmetric inhomogeneous 
C^l | strong magnetic field. However, there are some points which were not clearly explained. 

In fact, the arguments, which led us to the asymptotic formula, are based on a numerical 
study of the integral of Eq. (10), as we will see in the main part of this paper. The 
' aim of this work is to present this study in detail. 

o 

1 Introduction 

Many authors have dealt with the numerical study of the fermion-induced effective energy, 
in the presence of magnetic fields of the form of a flux tube, see Refs. [II El El IU El • An 
analytical study of the same problem in two dimensional Euclidean space can be found in 
Ref. 0. 

A similar topic, is the study of the gluon-induced effective energy in the presence of a 
colored magnetic flux tube jHlEI]. However, in this case there is an essential difference: the 
effective energy has nonzero imaginary part and so the system is unstable. 

It is worth mentioning the recently developed functional method in Refs. |U IE] , which 
is applicable even for magnetic fields which do not exhibit a special kind of symmetry. For 
example, this method can be applied to a system of two separated magnetic flux tubes. 

In Ref. j3] we presented an asymptotic formula (Eq. (49)) for the effective energy when 
the characteristic magnetic field strength B m = J B(x) ■ dS/ird 2 (or the magnetic flux $ of 
the tube) tends to infinity and the spatial size d of the magnetic flux tube is kept fixed. 
However, there are some points which were not clearly explained. In fact, the arguments, 
which led us to the asymptotic formula, are based on a numerical study of the integral of 
Eq. (10), as we will see in the main part of this paper. The aim of this work is to present 
this study in detail. The final asymptotic formula is given by Eq. (23). 
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In addition, in Sec. 4 we compared the asymptotic formula of Eq. (23) with the strong 
field limit of the derivative expansion and find agreement. 



2 Effective energy in the presence of a magnetic flux 
tube 

The 3+1 dimensional fermion-induced effective energy in the presence of a magnetic field is 
given by the equation: 

Ee fm+ i) = -^ T Tr\n(p 2 + m 2 f ) (1) 

where lp = (fi = 0,1,2,3) and — — ieA^. The gamma matrices satisfy the 

relationship {7^,7"} = 2g fiu and T is the total length of time. 

We have shown in Ref. j3] that a magnetic field with strength independent of z coordinate, 
directed towards z-axis, has renormalized effective energy per unit length Z^ 1 equal to 2 

<™U = I £(£(»> + m}) In(^J^) - I £ E [n} (2) 

where E{ n } are the eigenvalues of the planar operator (7 m .D m ) 2 (m = 1, 2), and {n} is a set 
of quantum numbers. Note that E{ n y > since the operator {^ m D m ) 2 is positive definite as 
the square of a hermitian operator. For the sake of simplicity we will drop the index (ren) 
in the rest of this paper. 

If the magnetic field is cylindrically symmetric with finite magnetic flux, the spectrum 
of the operator {p/ m D m ) 2 is continuous 3 . Thus we can set E{ n } = k 2 . 

In order to perform the summation over {n} that appears in Eq. (2), we need a density 
of states, given by the equation: 



1,3 \ / l,S 



1 dA(k) 

TT dk 



where J2i, s P^/l ee \k) is the density of states for free space, and 5^ s {k) is the phase shift which 
corresponds to I th partial wave with momentum k and spin s. The function A(k) is defined 
by the equation: 

A(*)= lim £ S ljS (k) (4) 

s,l=—L 

From Eqs. (2) and (3), if we drop the field independent term Eufts W aric ^ integrate by 
parts, we obtain 

1 r+°° ( k 2 + m 2 f \ 
E eff(3+1) = -^l kin [—^T 1 ) (A(fc) - c)dk (5) 



1 L Z is the length of the space box towards z direction. 

2 For the case of massless fermions (mf — 0) the renormalized effective energy is given by appendix A in 
Ref. 0. 

3 Besides continuous spectrum, zero modes may exist according to the Aharonov-Casher theorem |14j . 
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where c = lim^ +00 A(k). Our numerical study, in Ref. [3], shows that c = — 7t<f) 2 (<fi = 
e$/27r). This means that c is independent of the special form of the magnetic field we 
examine and depends only on the total magnetic flux of the field $. Also, in Ref. [3J, we see 
that the function A(k), when k — > +00, tends to c fairly rapidly, and thus, the integral over 
k in Eq. (5) is convergent. 

Eq. (5) is suitable for the numerical computation of the 3+1 dimensional fermion-induced 
effective energy in the presence of a magnetic flux tube. The calculation of the phase shifts 
is performed by solving an ordinary differential equation. For details see Refs. [SI E] • 



3 The asymptotic behaviour for the effective energy 
when B m — > +00 and d and mj are kept fixed 

The function A(k) depends on three quantities: k, B m = 2(p/d 2 and d. From these we 

construct the following two dimensionless quantities 4 : kB^I 2 and B m d 2 (=20). For dimen- 

_i 

sional reasons the function A(k) can be put into the form A(k) = G(kB m 2 , B m d 2 ). Setting 
_i ' _i 

A(k) = G(kB m 2 , B m d 2 ) and making the change of variable y = kB m 2 in Eq. (5) we obtain 

D "' '""y (HV 2 + m 2 f /B m ) - Hm)/B m )) (G(y, B m d 2 ) - c) dy (7) 



2tt 2 Jo 

In the strong magnetic field case (B m /m 2 » 1) we obtain: 

B 



£ e //(3+i) -»• - J o V^y [G{y,B m d )-c)dy 

r+oo . . 

_ 2^ ln (- Bm / m /) X y B ™ d2 ) ~ c ) d y 



There are many ways to achieve a large ratio B m /m 2 by changing the three independent 
variables B m , rnj and d. However, the most interesting case is when B m — > +oo and mj and 
d are kept fixed (see also Refs. PjJ El El)- This means that we keep the spatial size of the 
magnetic field configuration fixed and we increase the characteristic magnetic field strength 
B m or the magnetic flux <fi = B m d 2 /2. 

From Eq. (8) we see that the asymptotic behaviour of the effective energy when B m —>■ 
+oo (rrif and d are kept fixed) depends on the asymptotic behaviour of the integrals: 

h(B m d 2 ) = X y(G(y,B m d 2 )-c)dy (9) 

r+oo . . 

h{B m d 2 ) = ^ y\ny(G(y,B m d 2 )-c)dy (10) 
In the case of the integral of Eq. (9) it was possible to find an analytical expression: 

r+oo I r 

/ y(G(y, B m d 2 ) - c)dy = —— / d 2 xB 2 (x) (11) 



4 We have assumed the rescaling B m — * eB n 
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In order to give a proof for the above equation we will use Eq. (13) in Ref. [3] 

Tr(j m D m ) 2 = 5> {n} = ~l d 2 xB 2 (x) (12) 



{n} 

By using Eq. (3) for the density of states we obtain 



{n} \ 1,8 / 

If we drop the field independent term Yli s p\^ ee \k) and set again A(k) = G(kB m 2 , B m d 2 ) 
and y = kB^ 1 ^ 2 we find 



{n} 



o n dk 
= / fe(A(Ar) - c)rfA; 

7T JO 

9 /■+oo 

= — B m y{G{kB-! ,2 ,B m d 2 )-c)dy (14) 

7T JO 

Eq. (11) is a straightforward result of Eqs. (14) and (12). 
For convenience we write 

B(x) = B m B(x/d) (15) 

where the magnetic field B(x) has a characteristic magnetic field strength B m = 1 and range 
d = 1, Eq. (11) can then be written in the form 

j Q y(G(y, B m d 2 ) - c)dy = J d 2 xB 2 (x) (16) 

From the above relationship it is obvious that the integral J +o ° y(G(y, B m d 2 ) — c)dy is pro- 
portional to B m d 2 . 

I could not find an analytical expression, like that of Eq. (16), for the integral l2{B m d 2 ) 
of Eq. (10). However, a numerical study, which is presented in Fig. ^ suggests a linear 5 
asymptotic behaviour for this integral 6 , when B m — > +oo, independently of the particular 
shape of the magnetic field. Thus we can write 

r+oo , . 

J o y\ny(G{y, B m d 2 ) - c) dy -> cl 5 m d 2 (17) 

where the constant cl is independent of B m d 2 , and its value depends only on the special 
form of the magnetic field configuration. 



5 Our numerical study can not exclude the possibility of a weak logarithmic growing of l2(B m d 2 ) / B m d 2 , 
as we will see at the end of this section 

6 This implies that the asymptotic behaviour of the integral of Eq. (10) is the same with that of integral 

r+oc 



Jo °° y(G(y,B m d 2 ) -c)dy. 
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Figure 1: — (1/B m d 2 ) J °° yha.y(G(y, B m d 2 ) — c)dy as a function of B m d 2 for the magnetic 
fields B u B 2 , B 3 , B 4 . 



As we see in Fig. ^ we have performed numerical calculations for four magnetic field 
configurations with quite different kinds of inhomogeneity. 



B 1 (r) 

B 2 (r) 
B 3 {r) 
BJr) 



B r , 



1 



{r 2 /d 2 + l) 2 



B m exp(-— ) 

B m 6{d - r) 

B m 3(l- r -)6(d-r) 
d 



(19) 
(20) 
(21) 



We expect that our numerical work, presented in Fig. Q convinces the reader that the 
asymptotic behaviour of I 2 (B m d 2 ), as given by Eq. (17), does not depend on the kind of 
inhomogeneity of the magnetic flux tube. 

Thus, when B m — > +oo and mj and d are kept fixed, from Eqs. (8), (16) and (17) we 
obtain 



E, 



cl Bid 2 Bid 2 



e//(3+l) 



■ln(2WmJ) / d 2 SB 2 (x) 



We see that the logarithmic term dominates, so we can write 

Bid 2 



rpasympt 
^e//(3+l) 



24vr 2 



\n(B m /mi) / d 2 xB 2 (x 



(22) 



(23) 



5 



If, instead of the characteristic magnetic field strength B rn = &/ird 2 we use the parameter 
<p = e$/27r we see that the fermion-induced effective energy, for <p — > +00, is proportional 
to —(p 2 In </>. 

£:;s+\)~-^ 2 w (24) 

This asymptotic behaviour is correct, independently of the particular shape of the magnetic 
flux tube, as it was shown in the above numerical study (See Fig. |TJ). 

Note that the above asymptotic formula, of Eq. (23), is valid and for other cases of 
large ratio B m /m 2 . For example we note two: a) for rrif — > when B m and d are kept 
fixed and b) for d — > when rrif and <fi = B m d 2 /2 are kept fixed. For these two cases, the 
asymptotic formula, of Eq. (23), is obtained straightforwardly from Eqs. (8) and (16), and 
thus knowledge of the asymptotic behaviour of the integral I 2 (B m d 2 ) is not necessary. 

Finally, we note that our numerical study in Fig. ^is not sufficient to exclude a logarith- 
mic growing of l2{B m d 2 ) / B m d 2 . However this logarithmic dependence, even if it really exists, 
should be weak compared to the logarithmic dependence which appears in the asymptotic 
formula of Eq. (23), otherwise it would be clearly shown in Fig. ^and it is not. 

In order to observe a possible weak logarithmic growing we ought to perform computa- 
tions in an exponentially large range of values of B m . We could not execute computations for 
so wide a range, as the differential equation for the computation of the phase shifts (Eq. (69) 
in Ref. [3]) is stiff for large B m . However, we tried an estimate of this logarithmic dependence 
by fitting a curve of the form c\ + c 2 ln(B m d 2 ) to our data, as they are presented in Fig. [T] 
We find that c 2 = 0.0047 ± 0.0005 for B u c 2 = 0.0100 ± 0.0005 for B 2 , c 2 = 0.0062 ± 0.0005 
for B$ and c 2 = 0.0120 ± 0.0006 for B4. These results imply that, even if this logarithmic 
term c 2 \n(B m d 2 ) really exists, it would not contribute to the effective energy (as it is given 
by the asymptotic formula of Eq. (23) above) more than 5 per cent. 



4 Derivative expansion for B m —> +00 when d and rrif 
are kept fixed 

An approximate way for the computation of the fermion-induced effective energy is the 
derivative expansion ^21 EH] • This method gives accurate results for smooth magnetic field 
configurations. Note that the derivative expansion fails for the magnetic fields of Eqs. (20) 
and (21), as the magnetic field of Eq. (20) is discontinuous and the magnetic field of Eq. 
(21) has discontinuous first derivative. 

Our aim is to compare the strong field limit of the derivative expansion with the asymp- 
totic formula of Eq. (23). 

If we keep the first two terms of the derivative expansion, the effective energy, in the case 
of 3+1 dimensions, is given by the equation 

^//(li) = 4% w + 4+D t 5 ' ^ ( 25 ) 

where 
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For B m — > +00 from Eqs. (26) and (27) we obtain 

4% ~* -^HBJm)) J<PxB 2 (x) (28) 

4% - BJ.^n<'(-D J d^^g- (29) 

Eq. (29) can be derived straight from Eq. (27) if we assume that B m » m 2 (for details see 
Ref. |13j). In order to derive Eq. (28) from Eq. (26) we have used a similar way with that 
presented in Ref. [To] . 

Thus Eq. (28), which gives the strong field limit of the derivative expansion, is in agree- 
ment with the asymptotic formula of Eq. (23). 

5 Conclusions 

We emphasize that the arguments, which led us to the asymptotic formula of Eq. (23), 
are based on a numerical study. In particular, we computed numerically l2(B m d 2 ) for four 
magnetic field configurations with quite different kinds of inhomogeneity (see Eqs. (18), 
(19), (20) and (21)). Our results in Fig. UJ suggest that / 2 (-5 m <i 2 ) is proportional to B m d 2 
(for B m d 2 » 1) independently of the specific form of the magnetic field configurations. The 
consequence of this result is that also the final formula, which states that ^g//(3+i) ~ ~~ ^ m ^ 
for <p — > +00, is valid independently of the specific form of the magnetic field. 

Finally we showed, comparing with the asymptotic formula of Eq. (23), that the leading 
term of the derivative expansion gives the same asymptotic behaviour for the fermion-induced 
effective energy in the strong magnetic field regime even in the case (like that of the magnetic 
fields of Eqs. (20) and (21)) when the derivative expansion is invalid, because the next to 
the leading term diverges. 
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